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ABSTRACT 

Recent cosmological observations indicate that the present universe is flat and dark 
energy dominated. In such a universe, the calculation of the luminosity distance, d^, 
involve repeated numerical calculations. In this paper, it is shown that a quite efhcient 
approximate analytical expression, having very small uncertainties, can be obtained for 
d^. The analytical calculation is shown to be exceedingly efficient, as compared to the 
traditional numerical methods and is potentially useful for Monte-Carlo simulations 
involving luminosity distances. 
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1 INTRODUCTION 

The most recent cosmological observations indicate 
that t he present universe is fiat and vacuum domi- 
nated l|Komatsu et al.ll2009l ). In such a vacuum dominated 
space-time, the distance analysis requires computer inten- 
sive numerical calculations. Even though, computers to- 
day are very fast, efficient analytical calculation of distance 
, scales would be very useful for various types of Monte Carlo 
simulations. 

The most fundamental distance scale i n the uni verse is 
■ the luminosity distance, defined by (11 = \/ I//(47r/), where 
' / is the observed fiux of an astronomical object and L is 
its luminosity. Current astronomical observations indicate 
that the present density parameter of the universe satisfy 
S7a + = 1 with SIa ~ 0.7. Here f^A is the contribution 
from the vacuum and Qm is the contribution from all other 
fields. The distance calculations in such a vacuum dominated 
universe i nvolve repeated n umerical calculations and elliptic 
functions (lEisensteinl 1 1997l V 

In order to simplify the numerical calculations, IPenI 
l|l999h (hereafter Pen99) has developed quite an efficient 
analytical recipe. In this paper, we show another analytical 
method, similar in many respect to that of Pen99, that can 
be used to calculate the distances in a vacuum dominated 
fiat universe. 

Our analytical calculation is shown to run faster than 
that of Pen99 and has smaller error variations with respect 
to redshift (z) and f^A- 

Our recipe for calculating the luminosity distance is the 
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following {Hq is the present Hubble constant and c is the 
speed of light): 



dr — 



1 + 2 



[*(2:(0,nA))-*(x(2,nA))], (1) 

(2) 
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X = x{z, Qa) = ln(a -I- y/a"^ — 1), 



1 - f^A (l-f z) 



(3) 
(4) 



2 APPROXIMATION 

We first begin by analyzing how the scale factor, a{t) varies 
as a function of time t in a flat universe in which Qa 0. 
In this case, a{t) is given by l|Weinberd l200"il 

a = H^VLAa +HoQm-^, (5) 

where oq is the present value of the scale factor. The above 
equation is then immediately integrated into 

' cosh{3HotV^) - 1 • (6) 



1 ^^m 

2 



The scale factor is directly related to the z as, 
a _ 1 
ao 1 + 2 



(7) 



Let us define x = SJ/o^V^^a and indicate its present value 
by 2:0 = x{0, Q,a)- Then, equations [6] and [7] give 
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X — x{z,flA) ~ cosh ^ 



1 + 2 



Ha 



If we define a as follows 



a — a{z, Qa) = 1 + 2 ■ 



1 



l-^A (1 + ^)3 

and since a > 1 we can write x as 
X = x(z, Qa) = ln(a + •\/a^ — 1). 



(8) 



(9) 



(10) 



We note that a; is a monotonically decreasing function be- 
yond 1(0,0.7) = 2.42. 

We choose t he standard Robertson- Walker met- 
ric (|Weinberg||2008l) as the metric of the background space- 
time. With usual notation, this is 



ds 



1 — kr'' 



+ r'^{d9'^ + sin ed(f>^ 



(11) 



In the above space-time, we can use equation[S]to obtain 
r. A straightforward integration for a flat universe (fc = 0) 
yields, 



dx' 



[sinh ;|^]2/3 ■ 



(12) 



This integral can be evaluated in terms of hypergeomet- 
ric functions and related eUiptic integrals. But here we take 
an alternate, simple approach by defining a new function. 



*(a::) = lim 



dx' 



[sinh ^]2/3 ' 



(13) 



In the standard model the luminosity distance is defined 
as di = aor{l + z). Now we can use equation [T3]to write the 
luminosity distance as 

c 1 + z 



dL 



[*(xo) - *(x)] 



(14) 



Expanding ^I" in a series expansion to the 4th order, we 
find that 



*(x) = 3x''-'2 



l/3r,2/3 



X^ X^ 
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+ *(0), 



(15) 



where ^'(0) — —2.210. Now, equation 1141 reduces to the re- 
quired expression for the luminosity distance as 

l + z 



dL = 



[<f{xo) - *(s)] 



(16) 



3 ANALYSIS AND CONCLUSION 

In order to compare the method of Pen99 to ours, lets define 
the absolute relative percentage error as follows. 



AE : 



dr 



d-i 



X 100%. 



(17) 



Here d'^'i'"^ and dj,"" are luminosity distance values calcu- 
lated from approximate analytical methods and numerical 
method respectively. 

A comparison of AE for both analytical methods for 
SIa — 0.7 is shown in figure [T] Our method has a better 
absolute relative percentage error value for z < 1.0, 1.6 < 
z < 5.5 and z > 8.0 compared to that of Pen99. We note 
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Figure 1. The absolute relative percentage error (A_E) as a func- 
tion of the redshift for Qa = 0.7. 




Figure 2. Contour plot of absolute relative percentage error 
(AE) for the method of Pen99 with various z and Qa- 




Figure 3. Contour plot of absolute relative percentage error 
(AE) for our method with various z and Qa- 
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Figure 4. Histogram of running times of both Pen99 and our 
methods. 
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that the error in our method decreases steadily with redshift 
approaching < 0.014 % sd, z = 1100. In comparison, for high 
redshifts, Pen99 error always stays ~ 0.09% and does not 
decrease appreciably. 

A contour plot of AE based on the method of Pen99 
with various z and f^A is shown in figure O Relatively com- 
plicated distribution of variations in the AE can be seen for 
the parameter space characterized by z and Qa- However, 
a contour plot of AE for our method, which is shown in 
FigureO shows a smooth behavior over the same parameter 
space. 

In order to investigate the running time of the two ana- 
lytical methods we performed the following test. With z = 1 
and Qa = 0.7, we calculated the running time for 1 million 
calculations on a typical personal computer (Intel Core 2 
Processor, 2127 MHz, 1 GB RAM, IDlO Version 6.2 run- 
ning on Windows XP Service Pack 3) . Then we repeated the 
above process 100 times for the both methods. Histogram of 
both running time results are shown in figure[4l Our method 
is significantly faster than the method of Pen99. In addition, 
we performed the same test on the numerical method and 
found that our method is more than an order of magnitude 
faster. However, we note that the above test is hardware and 
compiler dependent and results may vary depending on the 
hardware and the compiler used. 

With less than 0.1% error, our analytical method be- 
comes quite desirable as the most interesting astronomical 
phenomena happen at z > 1 (JIa ~ 0.7). Furthermore, the 
analytical computation is more elegant and faster compared 
to traditional numerical computations invoked in connec- 
tion with calculations of distances in a vacuum dominated 
flat universe. 

Once we know the luminosity distance, it becomes a 
simple matter to evaluate the other distances such as the 
angular diameter distance or the proper distance. 
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